SMOOTHNESS OF DENSITY FOR THE AREA 
PROCESS OF FRACTIONAL BROWNIAN MOTION 



PATRICK DRISCOLL 

Abstract. We consider a process given by a two-dimensional 
fractional Brownian motion with Hurst parameter ^ < H < i, 
along with an associated Levy area, and prove the smoothness of 
a density for this process with respect to Lebesgue measure. 



1. Introduction 

Let B t := (Bj,B 2 ), t G [0, T] be two-dimensional fractional Brown- 
ian motion of Hurst parameter | < H < ^, and let 

{B m )t ■= {(B'Jt, (Bl) t ) 

denote the m-th dyadic approximation of B (as defined below in Sec- 
tion [22]). Define the area processes 

(1-1) (Am)t ■= \ 



(\B]) m d{B 2 s ) m - f (B 2 s ) m d{B]) 
o Jo 



and At := lim (A m ) t (where this convergence is almost sure - see 

m— >oo 

Theorem 2 of [7]). 

The main result of this paper is as follows: 

Theorem 1.1. Define the random process {Y}o<t<T, taking values in 
IR 3 , by 

Y = 0, 

Y t =(B t ,A t ). (te (0,7*1) 

Then for all t G [0, T], the density ofY t with respect to Lebesgue mea- 
sure is C 



The investigation of this process is motivated by the potential for 
fractional Brownian motion to be a useful driving signal in stochastic 
differential equations that model a wide variety of natural and financial 
phenomena; in particular, the presence of long-range persistence (for 
H > 1/2) or anti-persistence (for H < 1/2) makes fBm a natural 
candidate for a driving process in many scenarios. Several examples of 
such applications are included in [21] and [IT] . 

l 
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One area of interest in the study of stochastic differential equations 
is on finding sufficient conditions for existence and regularity of den- 
sities for solutions. More specifically, given some solution {Y t } to the 
equation 

d 

(1.2) dY t = y STx i (Y t )dZ t , 

i=l 

where {X{\ is some collection of vector fields and £ t is a Gaussian 
driving process on the space C([0,T],M. d ), it is natural ask whether £ 
admits a density with respect to Lebesgue measure. It is due to the 
celebrated theorem of Hormander (see, for example, Theorem 38.16 
of [20]) that, when our driving process is standard Brownian motion, 
our solution admits a smooth density so long as the set of vectors 
{X u [X^Xj], [[Xi,Xj],X k ] . . .} spans R d . 

In the case of fractional Brownian motion, one may no longer appeal 
to the types of martingale arguments used in proofs of the above result 
in the standard case. When H > -, the positive correlation of in- 
crements of sample paths results in better variational properties than 
those of Brownian motion, and so one may use Young's integration 
theory to attack the problem - existence of a density to a solution of 
(11. 2p under this condition is proven in [TJJ], and smoothness is proven 
in pp. When H < ~, one must turn to the rough path theory of T. 
Lyons (see [H]) in order to interpret (j 1.2ft in a meaningful manner. 
The connection between fractional Brownian motion with ^ < H < ~ 
and rough paths is investigated in [7] , and existence of a density in the 
case of | < H | is proven in [5]. As far as we are aware, Theorem 11.11 is 
the first positive result involving smoothness, and may give hope that 
similar results will hold in a more general setting. 

In order to prove Theorem 11.11 we appeal to the usual technique of 
Malliavin calculus. It follows from Theorem 5.1 of [15] that it is enough 
to show that Y satisfies the following two conditions 

1. Ye D°°; 

2. If 7 = DY(DY)* is the Malliavin covariance matrix associated 
to Y, then 

(det 7 )" 1 G L°°-(>V 2 ,P) := p|L J (W 2 ,P). 

We will begin by calculating the derivative DY explicitly, and from 
this Condition 1 will be proven in Proposition 13.91 Condition 2 will 
then follow by direct analysis of the Malliavin covariance matrix; see 
Proposition 13.131 and Corollary 13.151 
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2. Background 

2.1. Fractional Brownian Motion. A (one-dimensional) fractional 
Brownian motion {Bf 1 ; (t G [0,T])} of Hurst parameter H G [0, 1] is a 
continuous-time centered Gaussian process with covariance given by 

E(Bf Bf ) = R(s, t) := \ (> + - \t - sf H ) 

Our focus will be on fBm with 1/3 < H < 1/2; henceforth, we shall 
assume that such an H has been fixed and will drop the parameter from 
our notation whenever possible to do so without causing confusion. 

An n-dimensional fractional Brownian motion is a stochastic process 
{B t = (Bj, . . . , B™)\ t G [0, T]} is a continuous-time process comprised 
of n independent copies of one- dimensional fractional Brownian motion, 
all having the same Hurst parameter H. 

It is straightforward to check that the process B satisfies a self- 
similarity property; that is to say, the processes B at and a~ H B t are 
equal in distribution. By Kolmogorov's continuity criterion, the sample 
paths t H- B t are almost surely Holder continuous of order a, for any 
a < H - see Theorem 1.6.1 of [2] for details of this proof. 

2.2. Dyadic Approximation. For each m, we will let 

D m := {k2~ m T;k = 0,1,..., 2 m }. We define the m-th dyadic ap- 
proximator 7r m : C([0, T],R d ) — )• C([0,T],M d ) as the unique projection 
operator such that, for any given / G C([0, T], M d ), 

7r m f(t) = f(t), (teD m ) 

d 2 

^vr m /(t) = 0. (t(£D m ) 

In words, ir m f is nothing more than the piecewise linear path agreeing 
with / on the set D m . We will regularly use the shorthand notation 
fm '■= ffmf where convenient. Similarly, we will define the m-th dyadic 
approximation of fractional Brownian motion B m := ir m B; more ex- 
plicitly, 

(B m ) t := B t _ + (t- t.)2 m [B t+ - B t _], (0 < t < T) 

where t_ is the largest member of D m such that t_ < t and t + is the 
smallet member of D m such that t <t + . 

2.3. p- variation and Rough Paths. Let U be a Banach space (U will 
typically be R or M 2 ) with norm || • \\u, and P[0,T] denote the set of 
finite partitions of [0, T}. Suppose we are given a path / G C([0, T],U); 
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then for each 1 < p < oo and II = {0 = t±, t2, ■ ■ ■ , £jv = 1} £ V[0, T], 
one may define the quantities 

A*/ := f(ti) 
V P (f : H) := 




sup 

n€V[o,T] 

The norm || • || p is referred to as the p-variation norm; we shall define 
the space C P {U) :— {/ G C([0,T], U)\ \\f\\ p < oo}, which is a Banach 
space under || • It is easy to check that for given a and p such that 
a < -, any a-H61der continuous function is in C P {U). Also, one has 
that for any 1 < p < q, C p (U) C C q {U). 

Given / e C p (U),g G C g (U), where p and g are such that ^ + ^ > 1, 
one can develop the notion of integration of / against g in the following 
manner: if {U n := {ti}} C V[0, T] is a collection of partitions such that 
the mesh size |II n | tends to as n — > oo, we define 

,t #( n -) 



I f dg := lim V" /(q)A^ 

J v i=l 



where q G This limit is guaranteed to exist under the as- 

sumptions presented, and is independent of the choice we make of the 
family of partitions so long as their mesh size tends to zero - see The- 
orem 3.3.1 of [13] for further details. The element f dg is referred 
to as the Young's integral of f against g. This expression was originally 
formulated in [26] . We have the following estimate on the value of this 
expression (see Formula 10.9 of [26] ) 



(2.1; 



fdg-[f(0)-(g(T)-g(0))} 



<C\\f\\ p \\g\ 



where the constant C depends only on the values of p and q. 

Similarly, given some / G C([0,T] 2 , U), and partitions IT! = {s^}, 
n 2 = {tj} G P[0, T], we may define 

\jf ■= f(Si, tj) - f(s h tj-!) ~ /(Si_i, tj) + tj-!) 




V p (f : n 1; n 2 ) := V V ||Ay/||& 



( 2D ) : = sup y p (/ : n l5 n 2 

ni,n 2 eP[o,T] 
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As in the (one-dimensional) case above, we shall define the space 
Ci 2D) {U) := {/ G C{[0,T] 2 ,U); \\f\\ p 2D) < oo}. It is helpful to note 
that, for each / G C^ D \lI) such that /(0, •) = 0, then for each fixed 
8 G [0,T], /(*,.) G C P (U) and \\f{s,-)\\ p < \\f\\ p 2D \ since 

#(n) #(n) 

|A,/(a, -)| P = E l A ^( s ' •) " A ^(°' ■)!" 

i=l i=l 

#(n) 



<x;(iA i /( a) o-A i /(o ) .)r 

i Ai/(r, o - A*/( aj or) < (ll/ll^)" ■ 



i=l 



Trivially, one also has that for each / G C P (U), the function 
/<g) / : [0,T] 2 -»■ C/ defined by 

(J®/) (M) := /(*)/(*) 

is contained in D) (f/), and ||/ <g> f\\ ( p 2D) < \\f\\ 2 p . 

Given / G Cp 2Z ^(£7) and (7 G Cf D \lJ), where p and g are such that 
i+i > 1, the Young's integral off against g, denoted by f, Q r j 2 / dg, 
is defined by 

#(TI n ) #(* n ) 

/ f dg := lim V V /(q, dj) A^g, 

where, as before, {il n := {tj}}, := {sj}} C P[0,T] are collections 
of partitions such that the maximum mesh size |n n | V |\& n | tends to 
as n — > 00 and G 4/ £ ( s j-i> s j)- Existence of this limit 

under the given assumptions, independent of the choice of the family 
of partitions, is proven in Theorem 1.2 of [23], as is an estimate similar 
to that of the one-dimensional case: 

< C\\g\\W(\\f\\W + || /(0 , .)|| p + ||/(-,0)|| p + |/(0,0)|). 



/ dg 

[0,T]2 

It will be helpful to record here a pair of results relating the variation 
of paths with their linear approximations. 



Theorem 2.1 (Propositions 5.20 and 5.60 of [IT]). 

(1) Suppose x G C P (U), and let x m := ir m x be the dyadic approxi- 
mation to x as defined above. Then one has that 

|| •^m ||p — ^ || x||p. 
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(2) Suppose x G Cp D \lI), and let x m := 7T m x be the dyadic approx- 
imation to x as defined above. Then one has that 



I II (2D) < qp— In _| I (2D) 
I Jj m ||p _ II lip 



For further development of the theory, the interested reader may 
look in PS], and 



2.4. Gaussian measure spaces. Let (W, ||-||) denote a separable Ba- 
nach space. We will say that a measure P on W is Gaussian if there 
exists a symmetric bilinear form q : W* x W* — > M such that for all 
up e W*, 

J exp (vp(u)) dF(u) = exp (^~q{y, <P) 

Let B refer to the Borel cr-algebra on W; we will call the triple ( W, B, 
a Gaussian space. Define a continuous mapping J : L 2 (F) — > W by 



Jf := / uf(u) dF(u). 
Jw 

L 2 (P) 

Define % as the image of J restricted to the space W* ; this space 
may be equipped with inner product given by 

{Jf,Jg)n = (f>9)v*(F)- 

We will refer to Ji as the Cameron- Martin space associated to the 
Gaussian space (W, B, P). More information regarding the construc- 
tion of these spaces may be found in [3], [8], and [T2] . 

A well-known example of a Gaussian measure space is the one asso- 
ciated to Brownian motion where 

W := {w 6 C([0,T],R) : w(0) = 0} 

and the Gaussian measure P on W is the law of standard Brownian 
motion. Details of the construction of this measure may be found in 
and 

In this case, the Cameron-Martin space is 



U = Hi := yi E W : h(s) = J <f>{u) du; G L z [0, T] 
with inner product given by 

(h,k) Hl := f ' h'(s)k'(s)ds. 







A second example of a Gaussian measure space which is pertinent to 
the results described below, is as follows: let W be defined as above, and 
define the Gaussian measure P on W as the law of fractional Brownian 
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motion with Hurst parameter H < 1/2; then by following Proposition 
2.1.2 of [2], we have that % consists of functions of the form 
h(t) = f*K H (t,s)h(s) ds, where h G L 2 [0,T} and 



K H (t,s) := b H 



#-5 



(t-s) 



H~l) s^ H I (u- s)H - l -u H ~^ du 



where bu is some suitable normalization constant. The inner product 
of this space is given by 

(h,k) n := (h,k) L 2 [0tT] . 

For each fixed t G [0,T], the function R(t,-) = E[B t B] G U is a 
reproducing kernel for the space; that is to say, for any h G "H, we have 
the following: 

(h,R(t,-)) H = h(t). 

Let 5 refer to the space of cylinder functions; that is to say, random 
variables of the form 

FM^^H,..,^)), 

where / G C°°(R n ) with all partial derivatives having at most poly- 
nomial growth, and {fa, . . . , fa} C W*. Let D : 5 ->• S <g> "H* be the 
operator defined by the action by 

DF(u)k: = j t \t=oF (u + tk) 

n 

= 9if(fa(u>), fa(cu))fa{k) (ken). 

1=1 



For 1 < q < oo, we will let ID 1 ' 9 denote the closure of S with respect to 
the norm 

\\F\\ hq :=(E[\F\^}+E[\\DF\\U) 

One can naturally define an iterated derivate operator D k taking values 
in %® k ; from this we can define the seminorm 

||F|| fciff := (E[|Fn + ^E[||^F||^ fc ] 

V 3=1 
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and we will denote by B> k,q the closure of S with respect to || ■ j^g. Also, 
let 

D°° := p| p| B k ' q . 

fceNg>l 

Given some F G D 1,9 , we may define the Malliavin covariance matrix 
7 by 

7 := DF(DF)*. 

3. Proof of Theorem 11.11 
Again, we fix the Hurst parameter | < H < |. Let 

W 2 :=W8R 2 = WxM/ = {we C([0,T],M 2 ) : w(0) = 0}. 

On W 2 , one may construct a unique Gaussian measure P such that the 
coordinate process {B t } < t < T defined by 

B t (oj) = u(t) 

is a two-dimensional fractional Brownian motion with Hurst parameter 
H and P = Law(B). Our reproducing kernel Hilbert space in this 
instance is given as H 2 = "HcgilR 2 , where % is the Cameron-Martin space 
for one-dimensional fractional Brownian motion; for a general element 
h = (h 1 , h 2 ) e H 2 , the norm is given by ||/i||^ 2 = \\h l \\ 2 H + \\h 2 \\ 2 H . 
We shall fix 

< 31) pe (iT=W)' 

and define (following Section 5.3.3 of [TT]) the spaces 
W p := C°°([0,T\,R) n W IMIp ; 

W 2 := W p x w p ^ C oo ([0,T],M 2 )n>V 2lMlp . 

For reasons which will become apparent as we progress, it will be ben- 
eficial for us to declare that from here on out our process {B t } will 
be restricted to the probability space (W 2 , B^, P|w2); the details of 
this restriction are included in the Appendix. Most importantly, the 
Cameron-Martin space H 2 associated to the restriction of our measure 
is the same as the Cameron-Martin space associated to (W 2 ,S,P) as 
given previously. The following proposition shows that elements of % 
live within a smaller variational space, and will be used repeatedly in 
the sequel. 

Proposition 3.1. Let r :— Then the covariance kernel R associ- 
ated to fractional Brownian motion with Hurst parameter H has finite 
two-dimensional r-variation. Also, the associated reproducing kernel 
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Hilbert space H may be embedded in the space C r (M); furthermore, this 
embedding is a contraction. 

Remark 3.2. The assertions in Proposition \3. 1\ were originally made 
in p. 3363 of [5] or p. 13 of [TUJ. However, some technical gaps have 

been found in their proof that R G Ci 2D ^ ; as a result, we have included 
a direct proof of this statement within the Appendix (see Theorem \4-3\ 
below). 

Note that the above implies that "H 2 C Cp(M. 2 ), and, since r < p, 
H 2 C W 2 . Also note that Proposition 13.11 implies that R(s, •) G C r for 
each fixed s G [0, T]. 

Let Y m be the solution to the differential equation 

d(Y m ) s = X 1 (Y m ) s d(B 1 J s + X 2 (Y m ) s d (B 2 m ) s 

(^m)o = 

with vector fields given by X± := J| — \yj^,, X 2 := + Without 
a loss of generality, we may fix all of the processes in question at time 
T, and will elect to do so for the remainder. Explicit calculations show 
that 

(Y m ) T = ((B m )T, \ j\ B l)s d(B 2 m ) s - (B 2 J S d{Bl))j . 

As a result of Theorem 19 and Corollary 20 of [7], we have the 
existence of a process 

Yt := lim (Y m )T; (a.s. and inL 2 ) 

m— >oo 

we may denote this process suggestively as 

Y T = (b t ,± £ (Bl dB 2 - B 2 dBD^j = (B T , l -j\{B s ,dB s )^j 

(In this expression we have let a (a, b) := o}b 2 — a 2 ^ 1 .) 

3.1. Calculation of the Derivative of Y?. So that we might study 
the properties of the processes (Y m ) T and Y T , let us introduce the 
following expressions: given some u := bj x e\ + u 2 e2 G W 2 , we will let 
u be the element of W 2 defined by 



<jj := <jj 2 e\ — cu 1 e2 = Jcu, 
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where J is rotation by —ir/2. We will let q : "H 2 x "H 2 — > M denote the 
symmetric quadratic form given by 

(3.2) g(/i,fc) = lQf h(t) ■ dk(t) + k(t)-dh(t)^j 

= j\{t)-dk{t)+ l -(k{T).~h{T)). 



The integrals given in the definition of q are to be considered as Young's 
integrals; by our previous assumption that | < H < |, we have that 
2 = 4H > 1, and so q is well-defined. Furthermore, standard Young's 
integration bounds and Proposition 13.11 gives us that 

(3.3) |g(/i J A;)|<C r ||%||A:|| r <C'||/i||„||A;|| 7{ , 

and hence q is continuous in each variable. We may now rewrite our 
approximate solutions Y m in the following form: 

Y m = ^{B m ) T , -q(B m , B m )^j . 

As a consequence of the representation theorem of Riesz, we have the 
existence of a linear operator Q : H 2 — > H 2 such that (Qh, k) = q(h, k). 

Remark 3.3. Recall that r = j^; thus, our assumption on the value 
of p as made in ( \3.lh is sufficient to guarantee that the integrals above 
are well-defined, since this implies that - + - > 1. 

Proposition 3.4. Fix a G W p ; for each partition II = {ti}fL Q G 
V[0,T], define the vector Su G % in the following manner: 

N 

S n (-) ■= a ^ Wi> •) - ^fe-i' •)] > 
i=i 

where Ci G Then%— lim Su k exists, where {Yi^kLi C V[0,T] 

with | life | converging to zero as k — > oo; furthermore, this limit is 
independent of the family of partitions. We will denote this limit by 




a(t)R(dt,-). 



This limit satisfies the following properties: 
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(1) J a{t)R(dt. 

slant C > such that 

a(t)R(dt 



I[ot] 2 a ® a dR; hence, there exists a con- 



o 



<C\\a\\l\\R\\r\ 



(2) For each h G H, ( ^ a(t)R(dt, •), h\ = ^ a{t)dh{t). 

(3) (/ T a(t)^(dt,.)) (*) = Jo»(t)R(dt,s). 

Proof. First note that - + ^ > (1 — 2i7) + 2if = 1, which implies that 

(1) the Young's integral of a against R(-,s) for any s G [0,T] is 
well-defined, and 

(2) the 2D- Young's integral of a ® a against i? is well-defined. 
It follows from Theorem 2.1 of 123 that for each k, 



(3.4) 



1 1 Sri J 



#(n fc )#(n fe ) 

= Y Y a(c l )a(c j )(A i R(U,-),A j R(t j ,-)) H 
i=i i=i 

< Cllall^PH?^, 



where C is a constant depending only on p and r. Given any two 
partitions Il n = {si},U m = {tk} in the family, for c, G [sj_i, Sj], G 

||<Sn n - Sn m ||ft = l|Sn„||ft + ||Sn m ||ft - 2(Sn n , Sn m )ft 
#(n„)#(n„) 

i=l i=l 

#(n m )#(n m ) 

+ Y Y a(tZ fc )a(tZ ( )[Afcii2(tfc,ti)] 
fc=i j=i 

#(n n )#(n m ) 
- 2 ^ ^ «(ci)a(4) [Atf.i2(s 4 , t fc )] 

i=l fc=l 

n '^4°° / (a (8) a)(s,t) dR(s,t) 

J[0,T] 2 

+ / (a® a)(s,t) dR(s,t) 

J[0,T] 2 

-2 1 (a®a)(s,t)dR(s,t) = 0. 

J[0,T] 2 
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Hence, the completeness of % implies the existence of f Q a(t)R(dt, •). 
Since the 2D- Young's integral is independant of choice of partitions, 
one may also see from the calculation above that the limit of Su k is 
also independant of choice of partition, as claimed. Letting k tend to 
infinity in f!3.4[) and applying bounds as in Theorem 1.2 of [23] proves 
(1). For an arbitrary h G H, we note that 

/ / a(t) R(dt,-),h) = lim (S u ,h) n 
\Jo In l n l^° 

#(n) 

= lim y~] a(ci)(R(t i+1 , •) - R(t h -),h) 

1 1 i=i 
#(n) 

= lim a(ci) [h(t i+ i) - h(U)] 




a(t) dh(t), 



and so (2) holds. In particular, by setting h = R(s, •), (3) is a conse- 
quence of (2). □ 

Proposition 3.5. Let Q : H 2 — > H 2 be the bounded operator defined 
by 

q(h, k) = (Qh, k) u2 . 
Then the action of Q on elements of H 2 is given by 

(3.5) Qh:=-R(T,-)h(T)- [ h(t)R(dt,-). 

2 Jo 

Proof. Pick an arbitrary k G "H 2 . The inner product of k against each 
of the terms on the right hand side of (13.51) is given as 

(R (T, •) h (T) , k) u2 = h\T) (R (T, •) , k% + h 2 (T) (R (T, •) , k 2 ) H 
= h 1 (T)k 1 (T) + ~h 2 (T)k 2 (T) = ~h(T) ■ k(T), 
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and, as a result of Proposition 13.41 

/ h{t)R(dt,-),k) =([ J^i^Ridt,-),^ 

+ h 2 (t)R(dt,-),k 2 ^ 



h 1 (t)dk 1 (t) + [ h 2 (t)dk 2 (t) 
Jo 

h{t)-dk{t) = - [ h{t)-dk{t). 
Jo 



By combining these terms and comparing to (13.2)) . we see that the 
claim is proven. □ 



Proposition 3.6. Let Q : T-L 2 — > T-L 2 be the operator defined above. 

1. Q may be extended to an operator from W 2 into % 2 , which will 
also be denoted by Q; for any u G W 2 , 



Quo := \R{T, -)uj(T) - [ uj(t) R(dt, ■). 
* Jo 



2. Q is a bounded operator on W 2 . 



Proof. 1. That Q is well-defined as an operator on W 2 follows 
from Proposition 13.41 it is then of immediate consequence that 
Qu G U 2 for any u = (cu 1 , cu 2 ) G W 2 . 
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2. From this fact, along with the estimate given by ( 13. 4(1 . the value 
of || Qo; ||^2 for a fixed a; G W 2 may be calculated: 



IIQwII^, 



E 

i=l 
2 1 



-i?(T,-)^(T)- / Q*(t)R(dt, 



H 



i=l 



R(T,-)Co\T) 



H 



Co\t) R(dt, 



u 



(r(T,-)G)\T), £^{T) R(dt,- 



T 2H \Cj\T)\ 2 



E 

i=l 



T 



i=l 
+ 



cu J (T) / Q\t) R(dt,T) 
Jo 

T 2H \Co\T)\ 2 



4 



+ 



[0,T] 2 



(s,t) di2(s,t) 



a/(T)c^(t) i2(dt,T) 



< 



2 /T 2H \ 

Eii o< ie(V +||fl|| ^ )+||fl(r ' 0|| 



< 



2i3 



2|| J R||1 2D) 



which is finite by Proposition 13. 11 

□ 



Let us denote by QB the random variable taking values in R 2 : 



1 ' T 



QB:=-R(T,-)B T - / B t R(dt,-), 
* Jo 

where 5 := (B 2 , —B 1 ). 

We are now in a position to calculate the derivative of the process 
Yt- To begin with, for i = 1, 2 let us denote by i?£ the linear operator 
on R 2 with action given by 

R\h=(R{t,-)e h h)e i = ti(t)ei, 
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where {ei,e2} is the standard basis of M 2 . It is immediate that DB = 



Proposition 3.7. The process Yt has a derivative, DYt, taking values 
in the space of linear operators from H 2 into R 3 , with action given by 



DY T h = (R T h + R%h, (QB, /i)«a) a.s. 

Proof. Recall that Yt = {Bt, At), where At was defined as the almost 
sure limit of processes given by fll.il) . We have that Y T is continuously 
"H 2 -differentiable by [5J Proposition 3], and Corollaries 16 and 20 of [7] 
imply that E \{Y m )T — Yt\ 2 — > as m — > oo. We claim that DAtU = 
(QB, h); to prove this, it suffices to show that 



( \q(B m ,B r 



2 



Recall that the process B m was defined as the dyadic linear approxi- 
mator to our fractional Brownian motion B; similarly, we will denote 
by R m (u, v) the m-th dyadic approximation of the kernel R e % in the 
first variable; i.e., 



R m (u,v) := n m (R(-,v))(u). 
Since B{T) ■ B{T) = 0, it follows from (J3J| that 

q{B m ,B m )= / (B m )t ■ d(B m ) t 
Jo 



k=0 
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By the product rule (see Proposition 1.2.3 of [IE], for example), 



k=0 



i r 

Dq(B m , B m )h = -£ (HU) + • (B u - B h _ x 

r-T /-T 

h m (t) ■ d(B m ) t + (B m ) t -dh m (t) 



B t . + B t . 



h{T)-B T -2 / {B m ) t -dh m {t). 



T 



R(T,-)Bt-2 I {B m ) t Rm{dt,-),h 

/ K2 



Since = II^IIh f° r an y Hilbert space "H, we can use the above 

calculations to rewrite the left side of (13.61) as 



E 



( -R(T, .)Br - J (B m ) t R m (dt, 



E 



B t R{dt,-)- [ (B m ) t R m (dt, 



H 2 



Hence, to prove the claim is it required for us to show that 



(3.7) 



E 



B l t R(dt, 



(B l m ) t R m (dt, 



H 
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We note that for any h G H, 

(B m ) t R m (dt, •), h 



lim 

m— >oo 



H 



lim (i B rn)t k + (Bm)^) ( 

k=0 v 7 

lim Y](B m ) t ( h{t k ) - h(t k - 

m— >oo ' ' \ 

k=0 v 

+ (BmK., (h(t k ) - h(t k ^) 
2^ B t dh(t) = (J B t R(dt,-),h 



and so J^(B m ) t R m (dt, •) converges pointwise in % to J Q T B t R(dt, •). 

Applying Proposition 13.41 along with Theorem 12. 11 we may conclude 
that there exists a constant C for which 

2 



T 



B t R(dt, 



<C\\B\\ 2 p \\R\\i 2D \ 



T 



(B m ) t R m (dt, 



H 



<C\\Bf p \\R\\f D \ 



with the second inequality being independent of m. Hence, we may use 
Fernique's Theorem (Theorem 2.6 of [8]) to conclude that 

(3.8) I B t R(dt, ■)- I (B m ) t Rm(dt, •) < C\\B\\ 2 p < oo. 

JO JO H 

Hence, we may apply Dominated Convergence to conclude that (13.71) 
holds, as desired. 

□ 

Remark 3.8. In fact, we have shown something slightly stronger in 
the above proof. By changing the exponent on the left-hand side of 
(13.81) . we may conclude that \\D(Y m ) T — DY T \\^ n2 y , converges to zero 
in all LP ,j > 1. By applying the triangle inequality, we also find that 
\\DY T \\ m .eL°°-. 

Proposition 3.9. The random variable Yt is in 
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Proof. Corollaries 16 and 20 of [7] implies that E[|At| 2 ] < oo, and that 
A T = L? — lim (A m ) T . Hence, At is in the second-order Ito chaos; it 

follows from hypercontractivity (pp. 61-63 of jl8|, for example), that 
E[|Arp] < oo for all 1 < j • < oo. Combining this with the above 
remark, we find that 



|y T ||i . = E[\Y T \ j 



< E[\B T \ j ]+E[\A T \ j 



DY T\\\ n *y 



-E 



\DY T \\\ U 2y 



< oo. 



□ 



3.2. (det7) 1 G L°° (W 2 ,P). We begin by recording some more gen- 
eral results, which will be useful in proving integrability of (det7)~ J . 

Lemma 3.10. Suppose that X is a non-negative random variable such 
that, for each j > 1, there exists a constant Cj > for which 

E [e~ sX ] < C jS ~ j Vs > 1. 

Then X^ 1 G L°°~ . 

Proof. Fix j > 1. We note that for any k > 0, 

s^e'* 8 ds = k- j T(j), 

where T denotes the standard Gamma function. By letting k = X, we 
find that 



E[X~ 



1 



E 



ds 



1 



s^' -1 ! [e- sjr ] ds. 



Using the assumption given, we can see that this quantity is clearly 
finite. □ 

Theorem 3.11 (see Melcher [161 pp.26-27]). Let (W,B,P) be a Gauss- 
ian measure space with associated Cameron-Martin space %, and sup- 
pose $:WxW->lisa bounded non-negative quadratic form. Then 
the operator $ : H —> H given by 

<S>(h,k) = ($h,k) H 

is trace-class. In addition, if $ is not a finite rank operator, then 
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Proof. By Theorem 5.3.32 of [22], we have that for a set of independent, 
identically distributed standard normal random variables {Cn}^=i, the 
series B N := Yln=i ^nh n converges in W to B P-a.s. and in all LP,j > 1 
as iV — > oo, and 



Law Ij^Cnhn) =P. 



k n=l 



In particular, the fact that E||.B — _B ||yy — > implies that 

E\$(B N ,B N )\ ->E|$(B,B)|, 
and Fernique's theorem allows us to conclude that 



A 



N-roc 

n=l n=l 



J^(®h n ,h n ) = lim V$(/i n ,/i n 

n=l 

lim E [$(5^,5 



A"^ 

IL | ID , 

Af-s>oo 

E[$(£)] < CE[||£|p < oo. 



Thus, $ is trace-class. 

Suppose that <3> is not finite rank. Since $ is compact, there exists an 
orthonormal basis {h n ]'^L 1 C % for which = X n h n ; our assumption 
guarantees that #{n : A n > 0} = oo. Using this, it is easy to check 
that 



(B N ,B N ) = ($B N ,B N ) n = jr^i, 



n=l 



and so 



A 



*(B,B) = lim ]TA n ,e 

AT— >oo z — ' 



n=l 



We will let il^ := #{1 < n < N : A n > 0}; it is clear that K N 
Therefore, for each fixed N and positive s, 

( N 

E[exp(-s$ (B,B))] =E 



OO. 



< E 



exp is lim X n£. 

71=1 

N \ 
n=l / 



cxp 

N 



n 

71=1 



1 



2A n s + 1 



_ K N 

< Cms 2 



Applying Lemma 13.101 finishes the proof. 



□ 
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In order to apply Theorem 13. Ill we will explicitly calculate a formula 
for the determinant of the Malliavin covariance matrix associated to 
Y. 



Lemma 3.12. Given any a ^ 0, C 6 M mi „(R), and D e M n 
has that 

det a ^Tr ° D =a m (det(D-a- 1 C tr C)), 
where C tr is the transpose of C. 

Proof. This claim follows immediately when one writes 

I m a~ y C 
D — a~ l C tr C 



one 



dim 


c ' 











Qtr 


D 




Qtr 


In _ 





□ 



Proposition 3.13. Define the map 7 : W 2 M 3 (R) in the following 
manner: 

H L2 



7(w) := 



T 2H Io 



[(Qu(T))] tr \\Qu\\ 2 



(Qu) (T) 

H- 



Also, define the quadratic form $ on W 2 as follows: 

= T 4H \\Quf n2 - T 2H \Qu(T)\ 2 . 

Then 

(1) DYt(DYt)* = 7 a.s. 

(2) $ = det 7. 

Proof. 1. We begin by calculating the adjoint operator (Dy)*. 

Recall that i^/i = h^t)^. For a fixed A = (A 1 , A 2 ) 6 M 2 , 
/i G H 2 , and i = 1,2, 

<( J R^)*A,/ i > w2 = A-/i l (T)e 4 

= A^(T) = <A t J R(T,-)e J ,/i) w2 . 

Thus, one has that 

R T (PJ T )*\ = R^R(T, -)e 3 = 5 i:j T 2H \i. 
Suppose xeK 3 . 
((£>y T )*x, fc) = x • DYrk 

= x-(k(T),(QB,k)) 

= (x\x 2 )-k(T)+x 3 ((QB,k)) 

= x 2 ) + (R 2 T )*(x\x 2 ) + x 3 QB, k) 

= (x 1 R(T, -)ei + x 2 i?(T, -)e 2 + x 3 Q5, fc). 



SMOOTHNESS OF DENSITY FOR AREA OF FBM 



21 



Using this, we may now verify the claim: 



DY T (DY T )*x = DY T {x l R{T, -)e x + x 2 R(T, -)e 2 + x 3 QB) 



= DY T (x l R(T, -)ei) + DY T (x 2 R(T, -)e 2 ) + DY T (x i QB) 
= {T m x 1 , 0, (QB,x 1 R(T, -)ei)) 

+ (0,T 2H x 2 ,(QB,x 2 R(T,-)e 2 )) 

+ {x 3 QB l {T),x 3 QB 2 {T), (QB,x 3 QB)) 

T 2H I 2 QB(T) 1 
_(QB(T)f \\QB\\ 2 n2 \^ 



2. This follows immediately from Lemma [3.121 



Lemma 3.14. The quadratic form $ is positive semidefinite and has 
a trivial nullspace. 

Proof. The Cauchy-Schwarz inequality allows us to see that 



|(^)^(^)| 2 = |<^(^, 0, (Qc^)^^! 2 < ||^(^, OII^IKQ^)^!!, = ^^ik^)^ 



Hence, $ is non- negative and $ = if and only if Qu = cR(T, ■) for 
some constant vector c G 1R 2 . By the definition of Q, it would then 
follow that 



As a result of Lemma 31 of [9], one has that C^°(0,T) C H. Hence, 
we can conclude that (u)(t) — w(T)) + c is constant on [0, T], which 
implies that u(t) is constant as well. Thus u(t) = oj(Q) = for all 



Corollary 3.15. (<&) -1 e L°°-(W 2 ,P). 

Proof. Combining Lemma 13.141 and Theorem 13.111 gives us the desired 



□ 




which implies by Proposition 13.41 that 




te [0,T\. 



□ 



result. 



□ 



4. Appendix 



4.1. Two-Dimensional r- Variation for Covariance of Fractional 
Brownian Motion. Fix < H < 1/2 and T > 0. We will, as 
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usual, denote by R : [0, T} 2 — > R the covariance function for fractional 
Brownian motion with Hurst parameter H on [0,T]; that is, 

E[Bf Bf ] = R(s,t) := ~ [s 2H + t 2H -\t- s\ 2H ] . 

To this covariance function, we may associate a finitely addivite signed 
measure fin on the algebra generated by rectangles of the form 
{(a, b] x (c,d\ C (0,T] 2 } such that 

Hr((ci, b] x (c, d}) = R(b, d) — R(a, d) — R(b, c) + R(a, c) 

= -{\d- a\ 2H + \c- b\ 2H -\d- b\ 2H - \c - a\ 2H ) . 

It is easy to check the covariance of the process increments Bf — Bf 
and Bf — Bf is given by the //^-measure of the rectangle (a, b] x (c, d], 
or 

E[(Bf - Bf){Bf - Bf)] = m ((a, b] x (c, d}). 

Lemma 4.1. Suppose 0<a<b<c<d<T. Then the process 
increments Bf — Bf and Bf — Bf are negatively correlated; i.e., 
fi R ({a,b] x (c,d]) < 0. 

Proof. We begin by noting that the following relations hold: 

(1) {d - a) + (c - b) = {d - b) + (c - a). 

(2) (d - a) > [{d - b) V (c - a)] := max [(d - b), (c - a)]. 

(3) (c - 6) < [(d — 6) A (c — a)] := min [(d - 6), (c - a)]. 

Let f(a,(3) = a 2H + /3 2H . Then for any positive constant C, one can 
check that on the region {(a, j3) : a > 0, j3 > 0, a + (3 = C}, f has a 
maximum at f-) and decreases as either a or /3 are increased. Thus, 
/(d — a, c — 6) < /(d — 6, c — a), and the claim is proven. □ 

Throughout the sequel, we will make use of the fact that for H < ^ , 
one has the inequalities 

(x + y) 2H -x 2H <y 2H 

i i , . i 

for all x, y > 0. 

Lemma 4.2. Lei /i# 6e i/je measure defined above. Then for any in- 
tervals (a,b], (c, d] C [0,T], one /jas the bound 

M(a,b]x( C ,d])\<(b-a) 2H A(d-c) 2H . 
Proof. We will need to consider three possible cases: 
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(1) One interval is nested within the other, 

(2) the intervals partially overlap, or 

(3) the intervals are disjoint. 

Case 1. a < c < d < b. 

In this scenario, the claimed upper bound is clearly (d — c) 2H . Using 
this, we have that 

\fi R ((a, 6]xM)| = i|(d- a) 2H + {b- c) 2H -(b- d) 2H - (c - a) 2H \ 
<\{{d-af H -{c-af H ) 

+ 1 -((b-cr-(b- d r) 

<(d-c) 2H . 
Case 2. a < c < b < d. 
In this case, we know that 

MM] x MDI = ^\(d-a) 2H + (b-c) 2H -(d-b) 2H -(c-a) 2H \ 

<\{{d-a) 2H -{d-b) 2H ) 

+ \\{b-cf H -{c-af H \ 

<\{b-a) 2 »+ l -{{b-c) 2H V{c-a) 2H ) 
<{b- a) 2H . 
In a similar manner, 

\Hn((a, b]x(c,d[)\ = h(d- a) 2H + (b - c) 2H - (d - b) 2H - (c - a) 2H \ 
<\{{d-a) 2H -{c-a) 2H ) 

+ \\{b-cf H -{d-bf H \ 

<\{d-c) 2H + l -((b-crv{d-bD 



< (d-c) 2H . 
Case 3. a < b < c < d. 
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Here, we will use the concavity inequality twice to generate the de- 
sired bound. Firstly, we calculate that 



M(a, b] x (c, d})\ = \ | (d - a) 2H + (c - 6) 2H - (d - 6) 2 " - (c - a) 2 " | 

<\({d-af H -{c-af H ) 

+ \({d-bf H -{c-bf H ) 
< (d-c) 2H . 



In much the same manner, we find that 

|M(a, 6] x (c, = ^ | (d - a) 2H + (c - 6) 2 " - (d - 6) 2 " - (c - a) 2 " | 



<i((d-ar-(d-&r) 



+i( (c -ar-(c-6n 



<(6-a 



□ 



Theorem 4.3. Let r := i > 1. Then the function R has finite two- 
dimensional r -variation over [0,T] 2 ; more specifically, 

\\R\\ [2D) < (5T) 2H . 

Proof. Let 

n := {s := < si < . . . < s M ■= T}, 
:= {t := < ti < . . . < t N := T} 



be two partitions of [0, T\. Fix a j G {1, . . . , T}. We will let A be the 
unique integer such that sa-i < tj-i < sa, and L > A will denote the 
unique integer for which S£_i < tj < s L . 
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Figure 1. An example partition in s for a fixed strip 
tj—i <C t ^ tj. 



As usual, we define 

AijR := (A R ((si_i, Si] x (tj-.!, tj]) . 

Then 

M A-l 

(4.1) ^|A^r<^|A^r+|A^r 

i=l i=l 

i-1 M 

+ |A^r+|A^r+ ^ \A t3 R\ r . 

i=A+l i=L+l 

It follows from Lemma 14.21 that 

(4.2) lA^r < (tj - tj^), \A Lj R\ r < (tj - tj-t). 

Lemma [4.21 also implies that \AijR\ r < (sj — Sj_i); hence, we may use 
telescoping to bound the third term: 

L-l L-l 

(4.3) \ A v R \ T ^ Yl ( Si _ s *-i) = _ ^ ^ & _ 

i=A+l i=A+l 
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Let us now focus on the first and last terms of Equation ( 14. ip . Note 
that on each of these sums, Lemma 14.11 implies that A^i? < 0. We 
may use this fact along with Lemma 14.21 to see that 

.4-1 M 



i=l 



--L+1 
>A-\ 



M 



i J2 + E 1^1 

.i=l / \i=L+l 



E ^ 



2=1 



E A ^ 



i=L+l 



(4.4) 



= |/M(o,sa~i] x (t i _i,t j ])r 

+ \/i R ((s L ,T] x 
< 2(f i -f 3 -_i). 



Combining Equations (I4.ip - (l4.4p allows us to conclude that 



M 



i=l 

Hence, 

TV M N 

E i A v R \ r ^ E 5 fe - = 5T - 
j=i i=i j=i 

This completes the proof, since the two-dimensional r- variation of R is 
given as 



\R\e D) = I sup vviA^r <(5Tr. 

n,*67»[o,n n J 



□ 



4.2. Restriction of Gaussian Measures. At first blush, it may seem 
natural to have our process B have the classical Wiener space W 2 := 
C([0,T],R 2 ) as its sample space. However, doing so is not ideal, since 
many of the operators we will be considering are only defined on smaller 
spaces, such as the p-variation spaces. 

We begin with a general result regarding a-algebras. 

Lemma 4.4. Let X be any real separable Banach space and L be any 
Then \\-\\ x is a (£) - measurable if and only if Bx = o (£) . 
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Proof. It it easy to see that, in any case, cr(C) C Bx- Also, since ||-|| x 
is continuous it is always Borel measurable; therefore, if Bx = o~ (£) 
then \\-\\ x is clearly a (£) - measurable. 

Suppose that is a (£) - measurable; then for each x G cr(£), 
|| • — xo\\ x is also cr (£) - measurable, and x — > x — xo is a (£) /a (£) 
- measurable. From this observation, it follows that a (£) contains all 
balls in X. Since X is separable, every open subset of X may be written 
as a countable union of open balls. It follows, then, that a (£) contains 
all open subsets of X and therefore that Bx C o~(£). □ 

Theorem 4.5. Suppose (X, B = Bx, /•*) is a Gaussian probability space, 
and X is a linear subspace of X . Also let \\ ■ \\ x is a norm on X such 
that 

(1) The space (X , || • \\ x ) is a separable Banach space, 

(2) The embedding of X into X is continuous, 

(3) X G B and fi(X) = 1, 

(4) B := B x = {Af]X : A G B}. 

Then jx := fi\ x is a Gaussian measure and (X, B, ft) is a Gaussian prob- 
ability space. Furthermore, (X, /x) and (X, jx) share the same Cameron- 
Martin space H. 

Proof. Let R n /4 : X x X — > X x X is the rotation map defined by 
Rn/i(x,y) = (^(x-y),^(x + y) 



then by the rotational invariance of Gaussian measures (see, for exam- 
ple, Theorem 3.1.1 of [1]), proving the statement that jx is Gaussian is 
equivalent to proving that 



f(x, y) djx(x)djx{y) = fo R^/^x, y) djx(x)dfx{y) 

XxX JXxX 

for any bounded B x ^-measurable function /. Let / be such a function; 
since X is of full fx measure, we may extend / to an B x immeasurable 
function (which we shall also refer to as /) such that J XxX fdfxdfi = 
Ixxx fdfidfx (this extension may be done by setting a function equal 
to / on X x X and equal to zero on the complement, for example). 
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Then it follows that 




f(x,y) dfi(x)dfi(y) = / f(x,y) dfj,(x)dfi(y) 



f(x,y) dfi(x)dfi(y) 



f o R n/i (x,y) dfj,(x)dfj,(y) 



f o R n/i (x, y) dfi(x)dfi(y) 



f o R^ /4 (x,y) dp,(x)dp,(y) 



This proves the first assertion. 

To see the equivalence of Cameron-Martin spaces, we recall that 
J : L 2 (X,/i) X, defined by 



maps onto "H. Again, by virtue of fi being fully supported on X, we 
may extend any element of L 2 (X, ft) to an element of L 2 (X, //); thus it 



Remark 4.6. An alternate proof of the equivalence of Cameron- Martin 
spaces may be found in Proposition 2.8 of [8] . 

Let us now focus on restricting the law of fractional Brownian motion 
with Hurst parameter l/3<iJ< 1/2 to a variational space. The 
standard Gaussian space on which fBm is realized is (W, B, P), where 
W = {ue C([0, T], R) : w(0) = 0} and P = Lzw(B H ). Pick < e « 1 
and fix p := 1/H + e. Let < £ < T denote the evaluation map on 
W; i.e., = x{t) for any x G W. Since 



it follows that || ■ ||w is a cr({<j) t : < £ < T})-measurable function, and 
by Lemma [4.4} it then follows that c({(pt '■ < £ < T}) = Syy. Recall 
that we have defined the p-variation norm on W by 




is easy to see that J(L 2 (X,/i)) = J(L 2 (X, //)) = "H, as desired. 



□ 



W = SUp t , 
0<t<T 
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Recall that we have defined the space 

W p = {xGC oo ([0 1 T] 1 1):x(0)=0} IHIp . 

By Corollary 5.35 and Proposition 5.38 of [9], this space is a separable 
Banach space under the p-variation norm and contains all g-variation 
paths starting at zero for any 1 < q < p. Note that for x G W p , 
Holder's inequality implies that for any t G [0, T], 

\x(t)\ = \x(t) -x(0)\ 

< \x(t) -x(0)\ + \x(T) -x(0)\ 

< 2^ (\ x (t) - x{0)\ p + \x{T) - x(0)| p )^ 

p-i 

< 2 p F pj 



from which it follows that ||o;||>v < \\x\\ p , and so the embedding of Wq 
into W is continuous. Observe that we may rewrite the p-variation 
norm as 

/(#n) 



ne-p[o,r] 



i=l 



Thus, || • ||p is o"({0t|w p : < t < T})-measurable, which implies that 
cr(£) = £>w p - Furthermore, by Theorem 5.33 of [9], we know that the 
space W p is equivalent to 



#(n) 



x G W p : lim sup >. \ X (U) — x(t 
s ^°uer[o,T]:\n\<s ~[ 



i-lj 



If we now define 



#(n) 

a„(x) := lim sup > \x(U) — x{ti-i) 

n — 5.nn _ . , i < J 



n ^°° n&v[o,T}nN:\n\<± i=1 



then it follows that a p is a cr({(j)t\w p '■ < t < T})-measurable function, 
and that 

Wp = W p n {a p = 0} G B w . 
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Additionally, we may now use Lemma 14.41 to conclude that 
Bw P = <j{{<t>t\w P :0<t<T}) 

= {A n W p : A e <r({<j> t \w ■ < t < T})} 
= {An W p : A e B w }. 



Finally, we note that since the paths t H- Bf are a.s. Holder continuous 
of order /3 := H (l + ^f} 1 < H, each such path has finite q- variation 
for q = i = jj + § . So by Corollary 5.35 of [9], F(W P ) > P(>%) = 1. 
Thus, we may appeal to Theorem 14.51 to conclude that (W p , By\> , P|w p ) 
is also a Gaussian probability space, and that the associated Cameron- 
Martin space T-L coincides with the usual Cameron-Martin space corre- 
sponding to P on W. 
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